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1.0 INTRODUCTION 


This report describes the design, fabrication, and performance charac- 
teristics of two digital data signal enhancement filters which are capable 
of being inserted between the Space Shuttle Navigation Sensor outputs and 
the guidance computer. Commonality of interfaces has been stressed so that 
the filters may be evaluated through operation with simulated sensors or 
with actual prototype sensor hardware. 

The filters will provide both a smoothed range and range rate output. 
Different conceptual approaches are utilized for each filter. The first 
filter is based on a combination low pass nonrecursive filter and a cascaded 
simple average smoother for range and range rate, respectively. Filter number 
two is a tracking filter which is capable of following transient data of the 
type encountered during burn periods. A test simulator was also designed 
which generates typical shuttle navigation sensor data. Details of this 
design are included in the report. 

The design approach taken for both filters was to implement the filters 
as microprocessors. The first filter design was implemented as a special 
purpose device. The only component lacking in this filter to make it a true 
microprocessor is an arithmetic logic unit capable of add, subtract, logical 
operations, etc. This unit does contain an accumulator. 

The digital tracking filter was implemented as a true microprocessor. 

As a matter of interest, this unit could be programmed to perform the 
functions of the first filter as well as other filter concepts. 

As presently programmed, the first filter is capable of smoothing 
over 1, 2, 3, and 4 second intervals with sample rates of 16, 8, 4, and 
2 samples/sec, respectively. The output range rate word is computed based 
on the most recent range word and the seventeen previous range words. The 
smoothed range Word is based on the most recent range word and the fifteen 
previous range words. 


1-1 


The second filter is programmed to smooth over 2 , 4, 8, and 16 range 
words beginning with the most recent. Sample rates for any of the above 
options may be 2, 4, 8, and 16 samples/second. A smoothed range and range 
rate output are provided. Based on the computed range rate, both filters 
provide an update for the smoothed range word to correct for the time lag 
resulting from averaging over N words. 


2.0 FILTER ANALYSIS 


2.1 N-SAMPLE SMOOTHING FILTER 

The N -sample filter is a discrete, linear, time invariant, non- 
recursive, finite memory, polynomial smoothing filter. A review of the 
literature indicates that considerable work has been devoted to analysis 
of this smoothing technique. This analysis exists largely through the 
efforts of R. B. Blackman^ and J. D. Musa^v In this section, portions 
of this work are reviewed and applied to the problem at hand. 

Consider the random process R(nT) , where n is an integer and T is 
the sample period. The process can be thought of as comprising a signal 
component R(nT) and a noise component R(nT). It is assumed that R{nT) 
can be satisfactorily approximated by an rth degree polynomial in nT, 
R(nT). Also, assume that R(nT) is a random noise process that is wide 
sense stationary with respect to the sampling instants nT. This situa- 
tion occurs in the tracking of a moving object whose true position could 
be represented as an rth degree polynomial and whose measured position 
includes a random error. For significant filtering action to be possible, 
it is necessary that 

E[R(nT)] = 0 (l) 

and var [R(nT)J = var[R(nT)] = a R 2. That is, the signal must represent 
the d.c. component and noise is the time varying component over the 
time interval chosen for smoothing. 

Let $ R (iT) represent the autocorrelation function of R(nT) , where i 
is an integer. A discrete polynomial smoother of the pth order and rth 
degree is a filter that operates on R(nT) in such a fashion that the out- 
put C(nT) and the input R(nT) are related by 
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E[C(nT)] 2 = R ( P)(nT + t) 


( 2 ) 


for all n. The quantity t represents prediction time. If t is negative, 
the operation performed is interpolation and if t is zero the operation 
represents smoothing. The superscript (p) denotes the "pth derivative 
of the estimate with respect to nT." Thus for p equal zero, the output 
represents position data; for p equal one the output represents velocity 
information, etc. 

A linear, time- in variant smoother which is nonrecursive and has 
finite memory may be described in terms of the input-output relationship 


N-l 

C(nT) = ^2 W(1T) R[(n-i)T] 


(3) 


where W(iT) is the weighting function. Note that W(1T) is defined only 
at a finite number of points, that it is independent of the input, and 
that it is invariant with the time nT. No previous values of the output 
appear in (3), therefore the smoother is nonrecursive. The variance ratio 
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!d 

a R 2 


(4) 


is a measure of the efficiency of the smoother. Since R(nT) was assumed 
wide-sense stationary, C(nT) is also wide-sense stationary. This follows 
from the time invariance of the smoother. Therefore y 2 is not a function 
of time. 


If R(nT) has an autocorrelation function of arbitrary form, it may 
be shown using equations (1), (3), and (4) that 
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where W.. = W(iT) and s W(jT). If the power density spectrum of the 
noise is white, the autocorrelation function, <f>, assumes the form of 


* R [( 1 - j)T] 



and equation (5) becomes 


N-l 



1=0 


( 6 ) 


(7) 


However, the noise of interest here is not necessarily uncorrelated from 
sample to sample. Musa has shown that correlated noise may be represented 
by the wide-sense Markov process as a first order approximation, or by 
a linear combination of such processes as a better approximation. Either 
the autocorrelation function or the power density spectrum may be approxi- 
mated. 


2.1.1 Hide-Sense Markov Noise Model 

A wide-sense stationary, continuous random process will be called 
wide-sense Markov if it has the autocorrelation function 


$(t) ~ exp(-flx), t g 0. (8) 

A rigorous definition for wide-sense Markov noise may be found in Doob . 
The quantity Q, is called the "noise bandwidth." By using the eveness 
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property of autocorrelation functions, equation (8) may be written 
as 


®(t) = exp(-n|x| ) (9) 

For a discrete wide-sense Markov process with equally spaced samples, 
the autocorrelation function may be written as 


*(x) = exp(-fl| T | ) Cb T (T), (10) 

where Cby is the comb function defined by 

+« 

Cb T (x) = (11) 


The power density spectrum obtained by Fourier transforming (11) is 


S<f> = (2,f)2 + a 2 * T Cb l/T (f) (’2) 

where the * indicates convolution. 

Using the wide-sense Markov noise model reduces (5) to 

N-l N-l 

1,2 ' S X'i "j (13) 

i=0 j=0 

where 

ct = exp(-aT) . (14) 
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Applying this to the LM rendezvous radar, consider the system shown 
in Figure 1. 


White noise Sampler 


Radar 
Tracking 
Loop G(s) 

^ _ ff R 2 ^ 

Digital 

°C 2 

1 

Filter 



G(s) = - “ + 

s^ + 24.67s + 268.96 


Figure 1. LM Rendezvous Radar Tracking Loop Output With Smoothing Filter 


White noise is filtered by the radar such that the normalized power 
density spectrum at the input to the sampler becomes 

s(«) = 

where 

J (Gju) | ^dto = 12.10 (16) 

The closed loop transfer function is given by 



G( Ja) 


(15) 
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G(s) 


07 ) 


24(s + 11.2) < 

s 2 + 24.67s + 268.96 

Using Equations (15), (16), and (17), tne normalized noise spectral 
density at the radar output becomes 


s(w) 


26.58 


m 2 + 224.8 

oj 4 + 60.7u 2 + 72339 


(18) 


and is graphed in Figure 2 along with first and second approximations. 
The second approximation is given by 
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v=l 


(vnr + w* 


(19) 


where vn is the half power point in each approximation. 

An "optimum smoother" of the class being discussed has been investi- 
gated by Blackman. This specifies the filter whose weighting function 
yields the minimum possible value of y 2 . This work has been extended by 
Musa to include a linear combination of noise power inputs. The noise 
variance ratio of a linear discrete system for which the input noise 
power has been approximated may be expressed as 


* 2 • £\ 


1 + e 
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N - (N-2)e 
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( 20 ) 
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where N is the number of samples and T is the time between samples. Using 
this expression, the smoothing action of an N-sample averaging filter is 
graphed in Figure 3 for the LM rendezvous radar. The corresponding weight- 
ing function is given by 


N - (N-2)a ^ “ °» N_1 ) 

ir - ' iN-Zia ( f = '• 2 - — ■ N-2) 


( 21 ) 


and is shown in Figure 4. 

Referring now to Equation (2), for p = 1 . The filter output repre- 
sents velocity information. The variance ratio for this case is given 
by 


,2 _ 1 


J2TL. 


T (N-l) j [1 + n(N-l)] [2 + n(N-l)] + [1 - n 2 ]} 


( 9 2 ) 


where 


_ 1 - a 

n " ttt 


The corresponding weighting function i< 


(23) 
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(N-l ) j 

[ [l+n(N-l )][2+r] (N-l )] + 

[l-n 2 ]} 


' n 2 (N-l-2i) 


(N-l ) | 

J [l+n(N-l)][2+ n (N-D] + 
(l + ri)[l + n(N-2)l 

D-n 2 ]} 


(N-l) { [l+n(N-l )][2+n(N-l )] + [1-t, 2 ]} 


, (1 = 0 ) 

, (i = 1,2,— ,N-2) 
» (i = N-l) 


(24) 
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and both functions are plotted in Figures 5 and 6. 

2.1.2 Cascaded Simple Averages Smoothers 

The idea of obtaining a nonuniform weighting sequence of finite 
length by cascading uniform weighting sequences of shorter length origi- 
nated with J. W. Tukey 4 . The overall weighting sequence is the multiple 
discrete convolution of the cascaded simple weighting sequences. For 
example, the convolution of the simple weighting sequence 

12(l)j = jl, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 

with the weighting sequence 

j 6 0>] - jl> 1. 1. 1. 1. l| (26) 

is the weighting sequence 

12(1)1 * f(l) 


|l,2, 3,4, 5, 6, 6, 6, 6, 6, 6, 6, 5, 4, 3, 2, 1 1 . 


(27) 




This technique was extended by Blackman to include a differencing 
operation. This results in a simplification of the procedure for obtaining 
velocity or acceleration information. For example if 



then 



12 ( 1 ) 



|lj 1 , 1, 1,1, 1, 0, 0, 0, 0, 

0, 0, -1,-1, -1 , -1 , -1 , -1 | . (29) 
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This is equivalent to applying to the original time series a weighting 
sequence consisting of six consecutive ones, followed by six consecutive 
zeros, and followed in turn by six consecutive negative ones. 

2.1.3 S moothed First Differences 

Consider the following procedure. 

M Keep the most recent n T + 1 terms of the original time series 
in store. 

1.2 Generate each term of a second times series by taking the 
difference between the newest and oldest terms stored in 
1.1 and dividing by n-jT-. Keep the most recent n terms of 
this series in store (n 2 * n-j). 

1.3 Generate each term of the final output time series by summing 
the terms stored in 1.2 and dividing by n ? (see top half of 
Figure 7) . 

Thus, the second time series represents moving simple averages of n 1 con- 
secutive first differences and as such is a smoothed velocity. The whole 
scheme is equivalent to applying two moving simple averagings, in cascade, 
to the first differences of the original time series. In terms of 
Equation (29), the weighting sequence may be represented by 


^i n^n-gT j n 2^’ ^ n l“ n 2^ ’ n 2^"^| * (30) 
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These smoothers are an approximation of the optimum smoothers dis- 
cussed previously. The approximation for the first order smoother 
involves using only the values K, -K, and 0 for the weighting coefficients, 
where K is a constant. Referring to Figure 6, note that the optimum 
weighting function consists of a positive and negative triangle symmetrical 
about the middle of the time series. The approximation substitutes a 
rectangle of weight plus K for the first triangle followed by a weight 
of zero in the middle of the time series where the weights are small, 
and finally by a rectangle of weight minus K for the negative triangle. 

This technique reduces the amount of storage at the cost of a slight 
increase in p . 

The weighting functions of cascaded simple averaging smoothers for 
both the Oth and 1st order are as follows: 



l 

N 


and 


(31) 


(Osisf-1) 

w j " \ 0 (j * i * § N - 1) . (32) 

\ NT - (§"'N * i i N - 1) 

s 

Equation (32) is equivalent to Equation (30) where N is a multiple of 

3. The weighting function for the first order smoother is plotted in 
Figure 8. 
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Variance ratios for the Oth and 1st order smoother, are respectively : 


and 
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3N ( 1 -a) 
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- 2a 2N/ 3 - a N/ 3+ 2 ) 
[ n ( i - k )] 2 


(34) 


These variance ratios are plotted in Figures 9 and TO. 

The output noise variance ratio for white noise is given by 


2 

y 


27 

2 N 3 T 2 


0.14. 


(35) 


The first order smoother was investigated for providing a smoothed 
range rate. Comparing Figures 5 and 10 it is evident that the first 
order cascaded averages smoother provides smoothing action which is 
identical for all practical purposes. For ease of implementation, it 
this design for the range rate smoother was selected. 

It was desired to investigate the dynamic response of the range 
rate smoother under the shuttle functional requirement that range rate 

error be at most 0.3 fps. A computer program was written to evaluate 
the range rate error. Details of the program are contained in 
[5] - The simulations begin at a range of 10 nm with a closing rate 
of 30 fps. At 30 seconds into the r;n, a 20 second burn was initiated 
at a constant acceleration level of 0.5 fps. After the burn was complete, 
the runs were continued for an additional 20 seconds. The range rate 
error is plotted in Figure 11 for two sample rates. Note that the range 
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rate error is on the order of 0.17 fps during the burn period and zero 
for constant velocity. The abrupt change from zero to a constant error 
occurs because in the simulation the burn period start and finish coin- 
cided with the one second sample time of the guidance computer. Referring 
to Figure 10, for a two second smoothing interval, the reduction in range 
rate error by filtering action may be expressed as 


= yTr = 0.39. (38) 

That is, an input range rate error due to noise of 0.3 fps would be 
reduced to 0.12 fps. 

2.4 DIGITAL TRACKING FILTER 



The basic filter which was investigated for shuttle application 
is shown in Figure 12. 

In this filter, the delay time, T, is provided by a holding register 
which is clocked one sample time later to compute a new difference for 
the error. The previously smoothed range word is added to the smoothed 
error to compute a new range word. Multiplying the smoothed error by an 
appropriate constant results in a smoothed range rate. 

Two different designs for this filter were investigated in the 
initial analysis. In the first design, the smoothed error word consists 
of a running average which is used to update the output register each 
sample time. In the second design, the output register is updated after 
an N number of error words have been smoothed 

Performance of the loop insofar as stability and tracking are con- 
cerned may be investigated by obtaining the overall closed loop transfer 
function. The set of difference equations which relate the closed loop 
input and output in the time domain are: 
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FIGURE 12. DIGITAL TF 
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(39a) 






E(nT) = x(nT) - y(nT - T) 
and 

y(nT) = y(nT - T) + J-.[E(nT> + E(nT - T) + E(nT - 2T) 

+ . . . E[nT - (N - 1 )T] . (39b) 


Taking the z-transform of these equations, the closed loop transfer 
function may be shown to be 



1+z- 1 +z- 2 + .z~ (N-1) 

N - (N - Dz" 1 + z" 2 + . . ,Z N 


(40) 


In this configuration the loop filter has been assumed to be a 
straight n-sample averaging device whose output updates the output 
register each sample time. Stability analysis proceeds directly by 
determining the pole locations in this expression. If any poles fall 
on or outside the unit circle in the z-plane the system is unstable. 
Performing this analysis yields the information that the loop is unstable 
if N equals five or greater. For N equal four, the poles are located 
very near the unit circle. This analysis was verified independently 

with a computer simulation. Details of this program are contained in 
[5]. Figure 13 illustrates a simulated loop response with a 
ramp input. This simulation data has been verified analytically by 
inserting the appropriate transform for x(z) into Equation (40) and 
inverting the resulting equation for y(z). Note that for N equal five, 
oscillations about the mean never dec-ease because two poles are located 

outside the unit circle. 


v. .>• 
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The second version of this filter effectively employs a "skip 
sampler" inserted in the loop immediately preceeding the adder. This 
filter is illustrated in Figure 14, A correction factor is applied to 
-he smoothed range external to the feedback loop to compensate for lag 
associated with the averaging action in the filter. Similarly a correc- 
tion factor is applied external to the loop to insert the time base in 
the average error. This results in a smoothed range rate output. This 
filter was simulated and its range rate error is illustrated in 
Figure 15. Note that the range rate error is on the order of 0.5 fps. 

The large offset at the beginning of the run occurs because a large 
word was deliberately set into the smoothed range register at the start. 
These results, indicate that a digital tracking filter would be capable of 
producing high quality smoothed range and range rate data during vehicle 

acceleration. 

The tracking filter which was implemented is basically the configura- 
tion shown in Figure 14. Since this basic filter produces an output only 
after smoothing over N data words a multiplexing technique was devised 
which effectively implements N filters in parallel. With this technique 
each adjacent filter smoothes over N data words which Tags the previous 
set of N data words by one sample time. Thus, the filter presents a 
smoothed output word for each input word. 
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3.0 CASCADED AVERAGING FILTER DESIGN 
3.1 FIRST CANDIDATE FILTER 

The first candidate filter consists of a combination N-sample smoother 
which provides a smoothed range and a cascaded simple average smoother for 
providing range rate. A functional block diagram for the first filter is 
shown in Figure 16. During operation the smoothed range rate is computed 
immediately after the reception of a data word followed by the computation 
of the smoothed range word. The range rate computation implements Equation 
(32) and the smoothed range computation implements Equation (31). 

Operation of the filter may be understood by referring to the filter 
architecture illustrated in Figure 17. As indicated, the principal com- 
ponents are input-output registers, a 32 word memory, and an arithmetic 
logic unit. A programmed algorithm using ROM's is utilized to direct 
functional operations of the various logic units and a data bus provides 
I c * ata communications between the various units. 

The operational cycle begins when the filter raises the "data request" 
line. As soon as a "data ready" signal is received by the filter, the data 
request line is dropped and the data word entered into the input regester 
with subsequent storage into the memory. At this point the filter is pre- 
programmed to perform the following functional steps. The contents of the 
memory is scanned beginning with word #1 which is the most recent data word. 
This data word along with the previous five data words are summed together 
by the arithmetic logic unit which functions as an accumulator. Words 7 
through 12 are read out of memory but are not added by the accumulator so 
that these words are effectively skipped. Beginning with word 13 the next 
six words are read out of memory and sequentially subtracted from the 
accumulated sum which represents the first six words. In the event that 
subtraction through zero occurs, provision is made in the algorithm for 
complementing so that the resulting difference in the arithmetic logic 
unit is represented as a magnitude. At this point appropriate division is 
accomplished to enter the weighting function. This function is pre- 


3-1 


{ » 
























INPUT 

DATA 


Figure 17. First Candidati 



RANGE RANGE 

RATE 


Smoothing Filter Architecture 








determined by the setting of the switch position on the front panel. The 
computed range rate is stored in the output range rate register and the 
range rate data ready line is pulsed. At this time, the range rate is also 
stored in a scratch pad memory location for later use in updating the 
smoothed range word. 

The memory is connected as a carousel type such that after the eighteenth 
word, the pointer is automatically reset to word #1. A memory scan operation 
is begun and words #1 through #16 are summed by the arithmetic logic unit. 
Division by four is accomplished so that the average word is present in the 
accumulator. Since the smoothed word is delayed by eight sample times from 
the incoming data word, a correction factor equal to half the range rate is 
added or subtracted from the smoothed range word to update the computed 
range word. This is accomplished by retrieving the range rate word from 
the scratch pad memory location and either adding or subtracting depending 
on whether the range rate is increasing or decreasing. This word is stored 
in the output range register and the range data ready line is pulsed. This 
completes the operational cycle of the filter and it returns to the "idle" 
state until the next computation is initiated. A more detailed description 
of the filter's operation may be had by considering the algorithm which is 
programmed to operate the filter. 

3.1.1 First Candidate Filter Algorithm 

A flow diagram for the filter algorithm and its associated circuitry 
is shown in Figure 18 and 19, respectively. In Figure 19, the numbers in 
triangles serve to identify the ROM with the associated coding table. 

The coding which has been programmed in the ROM’s is contained in Tables I 
through IV. Numbers in circles refer to physical layout of the circuit 
board. 

Referring now to Figures 18 and 19, the filter remains in the "idle" 
state (step zero) until a "GO" signal is detected. The GO signal is derived 
from the filter's internal clock and has available repetition rates of 2, 

4, 8, and 16 Hz depending on the switch setting on the fron panle. This 
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Figure 19. First Candidate Filter Control Circuitry 










pulse sets the GO flip/flop (1/2 - 7474) so that Q equals a logic "1". 
Referring to Table I, ROM outputs Og, 0 fi , and 0 ? are coded for binary one 
whTch dials up position n on the 74151 multiplexer. The state of the GO 
F/F then appears at the »Y" output of the multiplexer. The presence of 
a logic ' 1" enables the 74163 counters so that the counters advance to step 
one on the next clock pulse. Alternatively, the presence of a zero at the 
"Y" output loads into the counters ROM outputs 0 Q through C, which, from 
Table I, are all zero's for step zero. 

After detection of the GO signal, the filter advances to state one and 
the data request line is raised. Referring to Table II, the data request 
ou-cput of ROM #2 becomes a logic "1" which is inverted so that the data 
request F/F is cleared, thus setting Q high. The filter then advances to 
step two and tests for data. In Table 1, ROM outputs Og, 0 g , and 0 7 are 
seen to be set to binary two. This dials up the third output of the multi- 
plexer which looks at the state of the data F/F. Upon receipt of a "data 
ready" from the simulator or navigation device, the data ready F/F is set 
such that Q equals logical "1". This result appears at the "Y" output and 
the ROM address counter advances the filter to state three on the next clock 
pulse and the remainder of the algorithm is executed. 

Rather than continue with a step by step description of the operation 
of the algorithm implementation, only specific comments will be made about 
the function of the remaining components of the circuit. 

Using the coding tables and the flow diagram any specific action of 
the circuit may be determined. The "A" counter is provided to enable the 
counter to perform the equivalent of a "DO LOOP". For example in step 
three, the 74157 which is essentially a two-pole switch routes the 0 r 
through 0 4 outputs of ROM #1 to the input of the "A" 'counters. These ROM 
outputs are coded for a binary nine which is leaded into the counter. 

Next, note that in step four the memory is read which places the contents 
of current data memory location on the data bus and an add command is given 
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to the arithmetic logic unit. These commands originate from ROM's #3 and 
4 (see Tables II and IV). At the completion of step four the filter 
advances to step five and the multiplexer dials up the carry output of the 
"A" counter. For C = 0, the address for step four is loaded into the two 
ROM address counters and the filter skips back to step four. This iteration 
will be performed six times at which time the "A” counter has reached count 
#15 which raises the carry out line and the filter advances to step #6. 

In Figure 19 note that one of the inputs to the multiplexer is the 
carry output of the arithmetic logic unit. When words 13 through 18 are 
being subtracted from the accumulated sum of words 1 through 6 subtraction 
through zero results in the complement of the difference word in the ALU. 

By "looking" at this carry during the subtract operation the filter detects 
this condition and complements the difference word so that the magnitude 
is contained in the ALU. At the same time the "SGN" F/F is set to denote 
that subtraction through zero has occurred. When the next word is read 
from memory the SGN = 1 conation directs the filter to simply add the 
magnitudes of the remaining six words. 


3 . 1.2 


First Candidate Filter Operation 

The Tirst candidate filter is capable of smoothing over 1, 2, 3, or 4 
second intervals. Corresponding sample rates for each of these smoothing 
intervals are 16, 8, 4, and 2 samples/sec respectively. The different 
smoothing intervals may be selected by the rotary switch on the from panel. 
A "DATA" light is included on the front panel and remains in the on state ’ 
while data is being transfered to the filter. In operation, after the 
filter is turned the correct smoothed words will be presented at the output 
after the first eighteen words. This number of data words are required to 
fill up the memory. For all subsequent operation, a correct smoothed data 
word will appear at the output for each input word. 


3* 1-3 First Candidate Filter Performance Data 

Verification of the operating characteristics of the filter was made 
using noisy simulated data. Figures 20 and 21 are data plots for the 8 
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LSB s of Input and output data. Noise was added to the 6 LSB's for both 
data runs. Since the same recorder was used for both input and output data, 

the data is not time synchronous but the smoothing properties of the filter 
are still apparent. 

In Figure 20 the data input consisted of noise only and four data 
runs were made using the four available sampling rates with corresponding 
sampling time. Using a true RMS voltmeter, the ratio of input to output 
noise variance was determined to be on the order of 10:1. The data run 
in Figure 21 is for a noisy ramp input. The simulator was set so that the 
range word increased three bits per sample time. 
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4,0 DIGITAL TRACKING FILTER DESIGN 
4.1 SECOND CANDIDATE FILTER 

The second candidate filter consists of a digital tracking filter which 
provides a smoothed range and range rate output. A functional block diagram 
for this filter is shown in Figure 14. As indicated previously, a multi- 
plexing technique was devised which permits N of these filters to be effec- 
tively implemented in parallel. The number N is controlled by switch 
settings on the front panel and may have values of 2, 4, 8, or 16. 

Operation of the filter may be understood by referring to the filter 
architecture illustrated in Figure 22. The principal components are input- 
output registers, two 32 word memories, and two arithmetic logic units. 

The memories are segmented into two blocks of sixteen words each. The 
memory units are identical with the exception that one memory board contains 
pull-up resistors for the data bus. 

The "DATA & ERROR" memory stores incoming data words in one block of 
16 words and error words are stored in the second block. In addition, this 
second block of memory is used as a "scratch pad.' 1 The range memory is 
used to store computed range words and is the unit which makes it possible 
to implement the N-parallel filter concept. 

Two arithmetic logic units (ALU's) are contained in the design. ALU #1 
is utilized strictly for performing addition, subtraction, and division. 

When computing error words, the results are stored in the error memory as 
a 15 bit error word with bit #16 (MSB) used as the word sign. When a nega- 
tive word results, the carry is monitored such that the word is automatically 
complemented so that the magnitude plus sign is stored. When other add or 
subtract operations are performed all sixteen bits are used as data. 

The second arithmetic logic unit is utilized is primarily as an accumu- 
lator which is also capable of left-shift and right-shift operations . Its 
primary use is to compute the average error word. The filter accomplishes 
this by scanning the error memory twice. On the first scan, ALU #2 sums up 
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an of the positive words. On the second pass, the negative words are 
sequentially subtracted from the accumulated sum. If subtraction through 
zero occurs, the contents of the accumulator are complemented and the re- 
maining negative words are added to the complemented word. A more detailed 
understanding of the filter operation may be had by examining the filter 
algori thm. 

4* 1 . 1 Second Candidate Filter Algorithm 

A flow diagram for the filter algorithm and its associated circuitry 
is shown in Figures 23 and 24 and the ROM codes are contained in Tables 5 
through 10. Referring to the flow diagram, the operation for data request 
and data acknowledge is identical to the first filter. In step #3 the 
current data word is loaded into the input register and also into the data 
memory and the filter advances to step #4. The previously computed range 
word is read from the range memory and loaded into the input register on 
the ALU #1 board. The output of this register forms the "A" input to the 
74181 's. In steps 5 through 7 an iteration is performed which scans the 
data memory beginning with the most recent data word and the difference 
between each data word and the range word is computed. These "error” 
words are stored in the error memory. When this operation is completed, 
the address pointer for the data memory will again be at the starting 
point which is the most recent data word. 

The address counter for the data memory is also used for the error 
memory. The "loop" starting with step #9 and ending with step #24 causes 
the error memory to be scanned. In this operation the word sign is examined 
and all of the positive error words are summed by ALU #2. The "loop" 
beginning with step #11 and ending with steps 28 or 30 subtracts the negative 
words from the accumulated sum. A complement operation is performed to 
facilitate subtraction through zero. At the conclusion of this operation, 
steps 13 or 14 perform a left-shift operation which divides the final sum 
in the accumulator to produce an error word which has been averaged over 
N words. This average error word has a value which is N times the range 
rate. 
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Figure 24* Second Candidate Filter Control Circuitry 
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In steps 17 and 18 this averaged error word is either added or sub- 
tracted from the current range word to produce the new range word. This 
word replaces the old range word in memory. Since the computed range word 
lags the actual range, a correction factor equal to N-l x R/2 is added or 
subtracted from the computed value of range. The remaining steps in the 
algorithm perform this operation. In this process, steps 31 and 32 per- 
form a divide by N operation which results in the actual range rate R. 

In step 38, the range memory is enabled so that the pointer advances 
to the next memory location. The data and error memory pointer is advanced 
during ~tep 40 or 41. These memory pointer advances effectively implement 
the N- parallel filter concept. 

4.1.2 Second Candidate Filter Operation 

The second candidate filter is capable of smoothing over 1,2, 3, 4, 
or 8 second intervals. Control of the smoothing interval is accomplished 
by varying either the sample rate and/or the number of words N to be used 
for data smoothing. The sample rate and number of words are determined by 
switch settings of four rotary switches on the front panel. Sample rate 
control is determined by the setting of the first switch on the left and 
its position is independent of the position of the remaining three switches 
The remaining three switches control the number of words N and must all be 
placed in the same position for correct filter operation,. 

In operation, the filter requires 2N input data words before it begins 
tracking. After this time, for a ramp input, the filter tracks exactly. 

For an abrupt step change in the input the filter will again require 2N 
input words to track the change. This characteristic is inherent in the 
algorithm. The filter is capable of processing error words of up to 15 
bits. Thus extremely large input data words can be tracked in a time frame 
corresponding to 2N input words. 


4-13 



Sec ond Candidate Filter Performance Data 

Verification of the operating characteristics of the filter was made 
using simulated noisy data. Figure 25 contains data plots for the 8 LSB's 
of input and output data. Noise was added to the 6 LSB’s for each data run 
As was indicated previously, the data plots are not time synchronous. 
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5.0 CONCLUSIONS AND RECOMMENDATIONS 


Based on recommendations contained in [5], two different filter types 
were constructed and evaluated. The results of this evaluation indicates 
that both filters will result in input/output noise variance reduction on 
the order of 10:1. In addition, the tracking filter demonstrated the 
ability to track large (up to a 15 bit word) discontinuities in the input 
data. The filter will reacquire lock on a data input change consisting of 
a full 16 bit word at a slightly greater expense in time. A computer aided 
evaluation will be necessary in order to obtain complete filter performance 
data. 

The implementation of both filters is based on extensive use of micro- 
programming techniques. Thus, the physical size, power requirements, and 
reliability of the filters may all be optimized by implementing the filters 
with commercially available microprocessor integrated circuits. The flow 
diagrams contained in this report may be implemented directly. 
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TEST SIGNAL GENERATOR 

A test signal generator was constructed to supply signals which can 
be used to evaluate the performance of the two filters. This signal generator 
supplies, under control of the filter under test, 16 bit range words which 
change at a rate determined by the setting of front panel switches. Pro- 
vision is also made for adding random noise to the output digital words. 

The overall operation of the test generator can be seen from an inspec- 
tion of the block diagram of Figure A-l. The range words are generated in 
an accumulator which increments the range word by a fixed amount upon each 
receipt of a data sample request from the filter under test. These range 
words are added to the noise words which are provided by an A/D converter 
which converts filtered analog noise from an external random noise generator. 
The resulting noisy range words are placed in an output register where 
they are available as 16 bit parallel words for use as input data to the 
filter under test. Two D/A converters are also included to convert both 
the output range words of the test generator and the output words of the 
filter to permit analog evaluations of filter performance. 

TEST GENERATOR OPERATION 

Schematic diagrams for each circuit board in the test generator are 
shown in Figure A-2 through Figure A-5 , as well as an overall interconnection 
plan as illustrated in Figure A-l. 

MAIN ALU BOARD 

The main ALU board provides the range word accumulator function. 

Four 74131 ALU's are connected to associated 74194 shift registers to form 
a 16 bit accumulator. The range word increment is set by the four switches 
S 2 through S g which are mounted on the front panel. The current range 
word appears on the output lines of the shift registers and is fed back 
as one input to the ALU's. The other ALU input is supplied by the range 
rate switches. When a clock pulse appears on the shift register clock line, 
the output of the ALU's is transferred to the range word output lines and 
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held until the next clock pulse. Initial clearing of the accumulator at 
power turn-on is provided by a one shot. Additional clearing of the 
accumulator is available from the front panel clear button which is con- 
nected so as to ground the clear line of all the shift registers? 

A/D CONVERTER BOARD 

This board contains an A/D converter for digitizing the output of the 
analog noise filter as well as facilities for adding the noise words to the 
output of the range word generator and for holding the resultant noisy data 
in an output register for use by the filter under test. 

When a data request is received from the filter under test, a "start 
conversion" signal is given to the A/D converter by a one shot which is 
triggered by the data request. This one shot also clocks the range word 
accumulator causing it to increment the range word. The six bit noise word 
output of the A/D converter is added to the current range word in a binary 
adder formed by four 7433^. The output of the adder is latched in the 
74194 output register by a clock pulse generated from the output of the 
"conversion complete" signal from the A/D converter. With appropriate 
delays, this same signal is used to clear the input register of the filter 
under test, to clear the data request flip/flop and to raise the data 
ready flag to the filter under test. 

LOW PASS FILTER BOARD 

An active low pass filter is used to filter the output of a random 
noise generate in order to simulate the low pass characteristic of the 
range tracking loop of a typical sensor. The active filter uses three 
operational amplifiers to develop a second order transfer function of the 
form 


H(s) 


S 2 + 2£w^S + 


(1) 
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The following sketch shows the configuration used. 



The transfer function is: 


V Q (S) 

W = S 2 + 


1 

1 

r 2 c 


s + 


Rp 


( 2 ) 


From Equation (1) it can be seen that the natural frequency, tu n , and the 
damping factor, £, can be controlled independently by appropriate changes 
in R-j and R,,. In particular, by comparing Equations (1) and (2): 

_ 1 

“n R 1 C 1 



( 3 ) 
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The values selected for C, R-j and R 2 in the test generator give co =62.5 
rad/sec., and 5 = 0.67. Other values are readily obtained by changing 
R-j and R^ appropriately. 

D/A CONVERTER BOARD 

It is often convenient to observe to the performance of the filter 
under test by means of an analog plot of the filtered and unfiltered range 
words. The D/A converter board contains two separate eight bit D/A con- 
verters for this purpose. Each converter is connected to the eight least 
significant bits of the range word being converted. Since the noise 
fluctuations are contained in the first six bits of the test words, con- 
version of the lower eight bits is sufficient to observe the behavior of 
the filter. If the generated range word is constantly increasing, the 
output of the D/A converter will reset itself after every change in range 
of 28 bits. The use of two converters permits simultaneous comparison of 
the filtered and unfiltered data. 
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Figure A-l. Test Signal 
Generator Block Diagram 
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APPENDIX B 

CASCADED AVERAGING FILTER 
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PIN NO. 


FUNCTION 


L&1&- 


ji 


1 Parallel Input Data -1 - MSB 


2 



-2 

3 



-3 

4 



-4 

5 


- 

-5 

6 



-6 

7 



-7 

8 



-8 

9 



-9 

TO 



-10 

11 



-11 

12 



-12 

13 



-13 

14 



-14 

15 

\ 

! 

-15 

16 

Parallel Input Data 

-16 

17 




18 




19 

Clock Out 



20 




21 

Shift Left Serial Input 


22 

Shift RT Serial Input 


23 

Data Request 


24 




25 

Data Ready 


26 




27 

Clear 



28 




29 

SO 



30 




31 

SI 



32 

6nd 



33 




34 





B-2 
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J2 


PIN NO. 
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7 
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10 
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14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 
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33 

34 


Parallel Range Rate Out 

A 


V 

Parallel Range Rate Out 
Paral lei Range Out 

A 


Parallel Range Out 
Range Ready 
Range Rate Ready 
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FUNCTION 
-1 - LSB 

A 


V 

-16 - MSB 
-1 - LSB 

* 


\K 

-16 - MSB 


t'Jf 
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appendix c 

digital tracking filter 
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Pin functions are the same as J1 and J2 for 
the cascaded averaging filter (see Appendix B) 
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